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ABSTRACT 


In  this  investigation,  a  continuum  mechanics  based  shear  deformable  shell  element  of  the 
absolute  nodal  coordinate  formulation  (ANCF)  is  generalized  to  a  laminated  shell  element  for 
application  to  the  modeling  of  fiber-reinforced  rubber  (FRR)  structure  of  the  physics-based 
ANCF  tire  model.  The  complex  deformation  coupling  exhibited  in  fiber-reinforced  composite 
materials  can  be  automatically  considered  in  the  shear  deformable  laminated  shell  element  using 
the  continuum  mechanics  approach,  and  the  element  lockings  are  systematically  eliminated  by 
the  assumed  natural  strain  and  enhanced  strain  approaches,  thereby  leading  to  a  locking-free 
shear  deformation  ANCF  laminated  shell  element.  Furthermore,  various  nonlinear  material 
models  can  be  considered  for  each  layer  in  a  way  same  as  solid  elements.  Using  the  ANCF 
laminated  shell  element  developed,  a  physics-based  ANCF  tire  model  is  developed  by 
considering  the  detailed  tire  geometry  and  material  properties.  The  experimental  validation  of  the 
tire  model  is  conducted  for  the  load-deflection  curve  to  ensure  that  the  fundamental  structural  tire 
properties  can  be  correctly  captured  in  the  ANCF  tire  model. 
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1.  INTRODUCTION 

An  accurate  modeling  of  the  complex  tire  geometry  and  the  anisotropic  material  properties  of 
tires  is  essential  to  the  tire  performance  evaluation  including  the  tire  contact  pressure  and  the 
braking/traction  and  cornering  forces.  Since  a  tire  consists  of  layers  of  plies  and  steel  belts 
embedded  in  rubber,  the  tire  structure  needs  to  be  modeled  by  cord-rubber  composite  materials 
and  various  fiber-reinforced  rubber  material  models  are  proposed  for  use  in  detailed  finite 
element  tire  models.  Since  Young’s  modulus  of  the  steel  cord  is  significantly  higher  than  that  of 
the  rubber  matrix,  mechanical  property  of  the  fiber-reinforced  rubber  (FRR)  is  highly  nonlinear 
[1].  In  particular,  the  tire  cross-section  property  that  includes  the  geometric  and  material 
properties  is  of  significant  importance  in  characterizing  the  normal  contact  pressure  distribution. 
Furthermore,  the  in-plane  shear  deformation  of  the  carcass  contributes  to  the  cornering 
characteristics  of  tires.  For  this  reason,  large-dimensional  high-fidelity  finite  element  tire  models 
that  account  for  the  tire  geometric  and  material  nonlinearities  are  developed  and  used  for  the  tire 
performance  evaluation.  However,  existing  finite  element  tire  models  cannot  be  integrated  into 
the  vehicle  dynamics  simulation  due  to  the  essential  difference  in  formulations  and  solution 
procedures  used  in  multibody  dynamics  and  nonlinear  finite  element  codes.  This  prevents  an 
integration  of  the  high-fidelity  tire  model  into  the  multibody  vehicle  dynamics  simulation  [2]  and, 
therefore,  the  structural  characteristics  of  tires  and  transient  tire  dynamics  are,  in  general, 
evaluated  using  different  computational  models  and  different  simulation  approaches.  To 
overcome  this  fundamental  and  essential  problem  in  the  tire  dynamics  simulation,  a  tire  model 
based  on  the  flexible  multibody  dynamics  approach  [2-4]  is  developed  using  the  absolute  nodal 
coordinate  formulation  (ANCF  [6,  7]).  The  in-plane  ANCF-LuGre  tire  model  developed  for  the 
transient  braking  analysis  allows  for  considering  the  nonlinear  coupling  between  the  dynamic 

I 

This  material  is  declared  a  work  of  the  U.S.  Government  and  is  not  subject  to 

copyright  protection  in  the  United  States.  Approved  for  public  release;  distribution  is  unlimited. 


UNCLASSIFIED:  Distribution  Statement  A.  Approved  for  public  release.  #26261 


Structural  deformation  of  the  tire  and  its  transient  tire  force  in  the  contact  patch  using  general 
multibody  dynamics  computer  algorithms  [4],  The  generalization  of  the  in-plane  ANCF  tire 
model  to  the  three-dimensional  model  requires  the  development  of  the  new  ANCF  shell  element 
suited  for  the  tire  model,  whieh  allows  for  modeling  the  nonlinear  fiber-reinforced  rubber 
materials  and  the  aecurate  three-dimensional  stresses  under  various  maneuvering  scenarios.  To 
this  end,  the  continuum  mechanics  based  shear  deformable  shell  element  of  the  absolute  nodal 
eoordinate  formulation  (ANCF)  [5]  is  generalized  to  a  laminated  shell  element  in  this  study  for 
applieation  to  the  modeling  of  fiber-reinforeed  rubber  (FRR)  strueture  of  the  physies-based 
ANCF  tire  model.  Furthermore,  a  physics-based  ANCF  tire  model  is  developed  using  the  shear 
deformable  laminated  shell  elements  such  that  the  high-fidelity  tire  model  can  be  integrated  into 
general  multibody  dynamics  computer  algorithms  for  ground  mobility  simulation. 

2.  CONTINUUM  MECHANICS  BASED  SHEAR  DEEORMABLE  ANCE  SHELL 
ELEMENT 


2.1  Kinematics  of  ANCF  Shell  Element 


As  shown  in  Fig.  1,  the  global  position  vector  r'  of  a  material  point  x'  =[x‘  y'  z'^  in  shell 
element  i  is  defined  as  [5] 


r'=r:(y,y)+z'A(v,y) 

CZ 


(1) 


where  r’^(x‘,y‘)  is  the  global  position  vector  in  the  middle  surface  and  dr‘(x‘ ,y‘)ldz‘  is  the 

transverse  gradient  veetor  used  to  deseribe  the  orientation  and  deformation  of  the  infinitesimal 
volume  in  the  element.  Using  the  bi-linear  polynomials,  the  position  vector  in  the  middle  surfaee 
and  the  transverse  gradient  vector  are  approximated  as  follows: 
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r:(V,y)  =  SL(x'-,y)e;, 


dz' 


(x'',y)=sux'',y)e; 


(2) 


where  =[5;i  5'I  5'I  5'l]  and 


i(l-^')(l  +  '7')  (3) 


where  =  2x‘  /  i"  and  rf‘  =  2y‘  /  w‘ .  £'  and  w'  are  lengths  along  the  element  x'  and  y'  axes, 
respectively.  In  Eq.  2,  the  vectors  and  e'^  represent  the  element  nodal  coordinates  associated 
with  the  global  position  vector  in  the  middle  surface  and  the  transverse  gradient  vector.  That  is, 
for  node  k  of  element  i,  one  has  =  r'*  and  e‘g  =  dr'^  / dz' . 

In  the  continuum  mechanics  approach,  the  elastic  forces  of  the  shell  element  are 
evaluated  as  a  continuum  volume  and  the  Green-Lagrange  strain  tensor  E  at  an  arbitrary 
material  point  in  element  i  is  defined  as  follows: 

E' =f((F‘'yF‘-l)  (4) 

where  F'  is  the  global  position  vector  gradient  tensor.  The  preceding  equation  can  be  expressed 
in  terms  of  the  covariant  strain  tensor  E'  as 

E''=(J')-^E'(J')-'  (5) 

where  J'  =dX‘/dx‘  and  X'  represents  the  global  position  vector  of  element  i  at  an  arbitrary 
reference  configuration.  The  covariant  strain  tensor  is  defined  as 

E'=^((j‘yj'-(j'yj‘)  (6) 

where  y=dr‘/dx‘  .  Using  Eq.  5,  the  strain  vector  e‘ =[£‘^  is 

defined  as 


3 

This  material  is  declared  a  work  of  the  U.S.  Government  and  is  not  subject  to 

copyright  protection  in  the  United  States.  Approved  for  public  release;  distribution  is  unlimited. 


UNCLASSIFIED:  Distribution  Statement  A.  Approved  for  public  release.  #26261 


£'=(Tr^£'  (7) 

where  e'  is  the  covariant  strain  vector  obtained  by  Eq.  6,  and  the  transformation  matrix  T'  is  as 
given  in  literature  [5]. 


2.2  Generalized  Elastic  Forces 

In  the  continuum  mechanics  based  shear  deformable  ANCF  shell  element,  element  lockings 
occur  due  to  the  use  of  low-order  polynomials,  resulting  in  erroneously  stiffer  bending  behavior. 
The  locking  in  the  continuum  mechanics  based  shear  deformable  ANCF  shell  element  includes 
the  transverse  shear  locking;  Poisson’s  thickness  locking;  curvature  thickness  locking;  in-plane 
shear  locking  [5].  These  lockings  are  systematically  alleviated  by  applying  the  assumed  natural 
strain  [8,  9]  and  the  enhanced  assumed  strain  approaches  [10,  11].  The  modified  strain  field  for 
the  continuum  mechanics-based  shear  deformable  ANCF  shell  element  can  then  be  defined  as 
follows  [5]: 


£  =  T  ^£  +  £ 


^EAS 


(8) 


where  the  covariant  transverse  normal  and  shear  strains  are  evaluated  by  the  assumed  strain 
approach,  while  the  other  covariant  strains  are  evaluated  as  compatible  strains  obtained  directly 
from  the  assumed  global  displacement  field.  This  leads  to  the  following  covariant  strain  vector: 

£  =  ls  s  f  (9) 

XX  ^yy  /  xy  ^  xz  ^  yz  _\  ^  ^ 


The  enhanced  assumed  strain  vector  =  \ 0  0l  in  Eq.  8  is 

|_  aa  yy  t  xy  zz  ^ 


defined  as  [10,  11] 


|Jo| 


To-"N(^)a 


(10) 


where  J(^)  and  Jq  are  the  global  position  vector  gradient  matrices  at  the  reference 
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configuration  evaluated  at  the  Gaussian  integration  point  ^  and  at  the  center  of  element  (^  =  0), 
respectively.  ^  is  a  vector  of  the  element  coordinates  in  the  parametric  domain  and  T,,  is  the 
eonstant  transformation  matrix  evaluated  at  the  center  of  element.  The  matrix  N(^)  defines 
polynomials  for  the  enhancement  of  the  strain  field  in  the  parametric  domain  and  a  is  a  vector  of 
internal  parameters  associated  with  the  interpolating  polynomials  N(§)  of  the  enhanced  strain 
field.  The  generalized  elastic  forces  of  the  shell  element  are  obtained  as  a  continuum  solid  using 
the  virtual  work  as  follows: 


Q'=-f  f-1 


dW‘(s‘) 

de‘ 


dV‘ 


(11) 


where  dV^  is  the  infinitesimal  volume  at  the  reference  configuration  of  element  i,  and  W  is  an 

elastie  energy  density  function.  The  eontinuum  meehanics  based  shear  deformable  ANCF  shell 
element  allows  for  considering  general  hyperelastieity  material  models  in  a  way  same  as  existing 
solid  elements. 


3.  LAMINATED  SHELL  FORMULATION 


3.1  Classical  Lamination  Theory 

In  fiber-reinforced  composite  materials  that  are  widely  used  in  many  engineering  applications, 
laminae  having  different  fiber  angles  are  bonded  together  to  produce  desired  material  properties. 
Since  many  laminae  are  stacked  at  different  fiber  angles,  the  complex  deformation  coupling 
between  the  extension,  shearing,  bending  and  twisting  occurs,  and  such  a  deformation  coupling 
characterizes  the  meehanical  behavior  of  fiber-reinforeed  composite  materials  [12].  In  the  first 
part  of  this  section,  the  macro-meehanieal  behavior  of  fiber-reinforeed  composite  materials  is 
overviewed  using  the  classical  lamination  theory. 
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In  the  Kirchhoff  plate  theory,  the  plane  stress  is  assumed  and  the  in-plane  strains  of  a  plate 
are  defined  as 


£p=<  +  ZK  (12) 

where  £°=[£’^  s^yy  ylyf  is  a  veetor  of  the  in-plane  strains  in  the  middle  plane,  and 
K  =  [/c’^  is  a  eurvature  veetor  associated  with  bending  and  twisting.  Using  a  linear 

orthotropic  constitutive  law,  the  in-plane  stress  vector  is  related  to  its  strain  vector  by  , 

and  the  material  moduli  is  given  as  [12] 

C^=R-'C^R-"  (13) 

In  the  preceding  equation,  the  transformation  matrix  R  is  a  function  of  the  fiber  angle  6  that 
defines  the  orientation  of  the  fiber  coordinate  system  o-12  with  respect  to  the  material  frame  o-xy 
of  the  plate  as  shown  in  Fig.  1.  This  matrix  is  defined  by 


cos^  0 

sin^  9 

sin  29 

sin^  9 

cos^  9 

-sin  29 

-sin  29 

sin  29 

cos26> 

(14) 


and  Cp  is  the  material  moduli  of  an  orthotropic  material  in  the  fiber  coordinate  system  as 


nil 


1122 


c 

c 


1122 


2222 


0 

0 


0 


0  c 


1212 


(15) 


where  =£■;(!- VJ2V21)  ^  =  V2i'E'i(l“U2V'2i)  ,  and  =0^^  • 

While  the  coupling  terms  between  the  normal  and  shear  strains  in  the  fiber  coordinate  system  are 
zero  as  observed  in  Eq.  15,  the  extension  and  shear  coupling  occurs  for  the  stress  and  strain  field 
defined  in  the  material  frame  and  the  coupling  terms  in  the  material  moduli  matrix  of  Eq.  13 
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are  not  zero.  That  is,  a  plate  subjected  to  a  uniaxial  load  produces  in-plane  shear  deformation  due 
to  the  extension  and  shear  coupling  [12].  Using  the  linear  constitutive  law  for  an  orthotropic 
material,  the  in-plane  stresses  of  a  lamina  are  defined  as 

a^=CX  +  zC^K  (16) 

The  force  and  moment  resultants  of  a  fiber-reinforced  composite  consisting  of  N  orthotropic 
laminae  can  then  be  defined  as 


N  = 


ZK.dz 


(17) 


and 


(18) 


where  N  =  [A^^  ,  M  =  [M^  ^  is  the  material  moduli  matrix  for  the 


^-th  layer,  and  z^  is  the  thickness  coordinate  at  the  upper  surface  of  the  k-th  layer.  It  is  important 

to  notice  here  that  the  force  and  moment  resultants  are  defined  as  forces  and  moments  per  unit 
length  [12].  Since  the  strain  and  curvature  are  not  function  of  the  thickness  coordinate  z  in  the 
Kirchhoff  plate  theory,  Eqs.  17  and  18  are  written  in  a  matrix  form  as  [12] 

(19) 


'N' 

A 

B' 

— 

P 

M 

B 

D 

K 

where  A,  =|;(C)„(Zi .  and  D,.  ='^(Cl\,(zl- zl,}l3  . 

k=l  k=l  k=l 

The  presence  of  the  matrix  B  implies  the  extension  and  bending/twisting  coupling  of  a  laminate. 

However,  if  each  lamina  above  the  mid-plane  is  identical  to  that  below  the  mid-plane  in  both 

geometry  and  material  properties  (i.e.,  the  laminate  is  mid-plane  symmetric),  the  matrix  B 

becomes  identically  zero  (i.e.,  5^  =0)  and  the  extension  and  bending/twisting  coupling  vanishes. 
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It  is  important  to  notice  here  that  the  extension  and  shear  coupling  still  exists.  Another  important 
case  is  a  balanced  laminate,  in  which  a  laminate  consists  of  a  pair  of  laminae  that  have  opposite 
fiber  angles  {+6  and  -6)  above  and  below  the  mid-plane,  regardless  of  the  stacking  sequence. 
This  eliminates  the  extension  and  shear  coupling.  That  is,  coupling  terms  of  in-plane  normal  and 
in-plane  shear  strains  in  matrix  A  are  identically  zero.  That  is,  for  a  two-layer  laminate  with  the 
same  layer  thickness  and  opposite  fiber  angle  {  +  1-9)  about  the  mid-plane,  the  extension  and 
shear  coupling  vanishes,  but  the  extension  and  twisting  coupling  in  B  matrix  exists.  This  causes 
twisting  deformation  of  a  laminated  plate  subjected  to  a  uniaxial  tensile  loading  [12].  The 
presence  of  such  a  complex  coupling  is  discussed  using  the  ANCF  laminated  shell  element  in  the 
last  part  of  this  Section. 

3.2  Generalization  to  Shear  Deformable  ANCF  Laminated  Shear  Element 


As  discussed  in  Section  2,  the  continuum  mechanics  based  ANCF  shell  element  is  formulated  as 
a  continuum  solid  that  accounts  for  the  three-dimensional  stress  state,  thus  the  complex 
deformation  coupling  exhibited  in  laminates  can  be  automatically  considered  without  special 
elastic  force  formulations.  That  is,  the  generalized  elastic  force  of  the  laminated  shell  element 
that  consists  of  N  layers  can  be  defined  as  follows: 


5£' 


ik  \ 


dV! 


(20) 


In  the  preceding  equation,  the  integration  interval  for  the  k-th  layer  in  the  thickness  direction  is 
from  to  Zi, .  In  other  words,  the  element  generalized  elastic  forces  are  evaluated  layer  by 

layer  and  the  resulting  generalized  elastic  forces  of  each  layer  are  simply  added  together  to 
define  the  elastic  force  vector  of  the  laminated  shell.  It  is  important  to  notice  here  that  there  is  no 
restriction  in  material  models  considered  in  each  lamia,  despite  the  fact  the  orthotropic  material 
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law  is  the  most  popular  material  model  used  for  reinforced  composite  materials.  Two  Gaussian 
integration  points  are  used  along  the  thickness  when  the  elastic  forces  of  each  layer  are  evaluated. 
Orthotropic  Saint- Venant-Kirchhoff  Material  Model  For  an  orthotropic  Saint- Venant- 
Kirchhoff  material,  the  material  moduli  =  d^W  /  dsyds^.,  of  an  orthotropic  lamina  in  the 
material  frame  is  defined  as  [13] 

=  (b''  -aJCb^'  -a.Xb^  -aJCb'  (21) 

where  (J')-'=[b^  b^  b'l  and  is  the  tangent  material  moduli  defined  using  9  material 
parameters  in  the  fiber  coordinate  system  [Uj  aj  as  shown  in  Fig.  1,  where  the  direction  of 


fiber  is  defined  along  the  coordinate  1 .  The  material  moduli  in  the  fiber  coordinate  system 


are  given  as  follows  [13]: 


[C'^“]  = 


^1111 

^1122 

0 

^1133 

0 

0 

^1122 
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^2233 

0 
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0 
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^2233 

0 
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0 

0 

0 

0 
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0 

^2323 
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0 

0 

0 

0 

0 

^1313 

(22) 


where 


C""  =£i(1-V23F32)/A 

C"''"=£2(1-V'i3V3i)/A 

^3333 

C"^^=£3(V2,+V3,V23)/A 

=  £3  (Vj3  +  Vj2V23)  /  A 

^2233 

y^i2i2 

C  —  C7j2 

y^2323 

^  “  ^23 

^1313 

and  A  =  1- V12V21  -  V23V32  - 

-2V21V32V13. 

(1  ^12^21 )  ^  ^ 

=  £'2(V32+Vj2V3j)/A  (23) 

=  Gi3 


Mooney-Rivlin  Material  Model  For  modeling  incompressible  materials  such  as  rubbers, 
Mooney-Rivlin  material  model  is  widely  used.  The  energy  density  function  is  defined  as  [14] 

W  =  q(l, -3) +  C,(T, -3) +  ^(J -If  (24) 
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where  Cj  and  are  material  constants,  /j  =  /j  /  ,  4  =  /2  /  and  J  =  ,  where  l\, 

h  and  h  are  invariants  of  right  Cauchy-Green  tensor  [14].  ^  is  a  bulk  modulus.  The  second 
Piola-Kirchhoff  stress  tensor  S  is  obtained  by  differentiating  the  energy  density  function  W  with 
respect  to  Green-Lagrange  strain  tensor  E  as 


dW 


S  =  ^  =  2C,(/3) 


-1/3 


r 


aE 


1 


'1-1 


I  —  LC 

V  3  '  y 


+  2C2(/3) 


-2/3 


/jl-c  —  +/f/(J-l)C“‘  (25) 

V  3  J 


Equations  of  Motion  Using  the  principle  of  virtual  work  in  dynamics,  the  equations  of  motion 
of  the  shear  deformable  laminated  shell  element  i  can  be  expressed  as 


M'e''=Q'(e'',a')  +  Ql(e'',e'',0 


(26) 


where  vectors  Q'^  and  are,  respectively,  the  element  elastic  and  external  force  vectors;  and 
the  matrix  M'  is  the  constant  element  mass  matrix  defined  by 

(27) 

k=\  » 


where  is  the  material  density  of  k-th  layer  at  the  reference  configuration.  The  internal 

parameters  a‘  introduced  for  the  enhanced  assumed  strains  are  determined  by  solving  the 
following  equations  [10,  11]: 


d£ 


EAS 


6a‘ 


dW‘(E‘) 
6£‘ 


dV’  =  0 


(28) 


It  is  important  to  notice  here  that  the  preceding  equations  can  be  solved  at  element  level  for  the 


unknown  internal  parameters  a'  using  the  procedure  presented  in  the  literature  [5]. 


3.3  Comparison  with  Analytical  Solutions 

To  validate  the  ANCF  laminated  shell  element  presented  in  this  Section,  a  uniaxial  tensile  test  of 
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a  two-layer  laminate  at  +1-9  fiber  angles  is  considered  as  a  benchmark  problem  as  shown  in 
Fig.  2  [12].  The  length,  width  and  thickness  of  the  laminated  plate  are  2.0  m,  1.0  m,  and  0.01  m, 
respectively.  The  thicknesses  of  both  layers  are  same  (i.e.,  0.005  m),  but  the  sign  of  the  fiber 
angles  of  the  upper  and  lower  layers  are  opposite  as  shown  in  Fig.  2.  The  tensile  distributed 
uniaxial  load  of  500  N/m  is  applied  in  the  X  direction.  As  discussed  in  Section  3.1,  the  extension 
and  bending/twisting  coupling  described  by  B  matrix  in  Eq.  19  are  not  zero  in  this  laminate, 
causing  the  warpage  (twisting)  under  the  uniaxial  tensile  loading  [12].  It  is  also  important  to 
notice  in  this  example  that  the  extension  and  shear  coupling  exists  in  each  lamina.  However, 
shear  deformation  of  the  upper  and  lower  layers  developed  by  the  uniaxial  tensile  load  are  same 
in  magnitude,  but  opposite  in  direction.  For  this  reason,  the  shear  deformations  of  both  layers  are 
canceled  out  and  no  in-plane  shearing  occurs  at  the  laminate  level.  That  is,  the  extension  and 
shear  coupling  term  of  A  matrix  of  the  laminate  in  Eq.  19  is  identically  zero  in  this  problem.  To 
demonstrate  this  fundamental  coupling  behavior  of  the  laminated  composite  material  using  the 
ANCE  laminated  shell  element,  twisting  angles  of  the  two-layer  laminated  plate  are  presented  in 
Fig.  3  as  a  function  of  the  fiber  angle  0 .  In  this  figure,  results  obtained  by  the  ANCF  laminated 
shell  element,  the  laminated  solid  shell  element  in  ANSYS  (SOLSH190),  and  the  analytical 
model  discussed  in  Section  3.1  are  compared.  In  the  analytical  model,  the  warpage  is  defined  by 
[12] 

^  =  (29) 

and  the  curvature  vector  k  =  is  determined  by  solving  the  following  equation: 


A  B 

'N' 

p 

— 

K 

B  D 

M 
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for  N  =  0  0]^  and  M  =  [0  0  0]^  .  It  is  observed  from  Fig.  3  that  the  twisting 

deformation  obtained  by  the  three  different  models  agree  well,  and  the  warpage  developed  by  the 
uniaxial  tensile  load  applied  to  the  two-layer  laminated  plate  is  well  predicted  by  the  ANCF 
laminated  shell  element  developed  in  this  study.  It  is  important  to  notice  here  that  the  sign  of  the 
twisting  angle  of  the  composite  plate  changes  at  the  fiber  angle  of  0  =  54.1  degree,  and  this 
important  fiber  angle  is  also  correctly  predicted  with  the  ANCF  laminated  shell  element. 

4.  APPLICATION  TO  PHYSICS-BASED  ANCF  TIRE  MODEL 

A  tire  has  a  complex  structure  that  consists  of  layers  of  plies  and  steel  belts  that  are  embedded  in 
rubber,  thus  an  accurate  modeling  of  the  complex  tire  geometry  and  the  anisotropic  material 
properties  is  essential  to  the  tire  performance  evaluation  including  the  tire  contact  pressure  and 
the  braking/traction  and  cornering  forces.  While  the  in-plane  tire  belt  deformation  can  be 
modeled  by  an  equivalent  material  model  [3,4],  such  a  simplified  material  model  cannot  be  used 
for  predicting  the  overall  tire  structural  deformation  in  the  three-dimensional  analysis.  This  is 
attributed  to  the  fact  that  the  tire  section  property  in  both  geometry  and  material  is  of  crucial 
importance  in  characterizing  the  contact  pressure  distribution.  This  necessitates  various  fiber- 
reinforced  rubber  material  models  that  can  be  integrated  into  high-fidelity  finite  element  tire 
models  [15,16].  However,  despite  the  fact  that  accurate  solutions  can  be  obtained  using  existing 
finite  element  tire  models,  difficulties  arise  when  they  are  integrated  into  the  vehicle  dynamics 
simulation  due  to  the  essential  difference  in  formulations  and  solution  procedures  used  in 
multibody  dynamics  and  nonlinear  finite  element  codes.  This  prevents  an  integration  of  the  high- 
fidelity  tire  model  into  the  multibody  vehicle  dynamics  simulation.  That  is,  the  structural 
characteristics  of  tires  and  transient  tire  forces  are,  in  many  cases,  evaluated  using  different 
computational  models  and  different  simulation  approaches.  For  this  reason,  a  physics-based 
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absolute  nodal  coordinate  formulation  tire  model  using  the  shear  deformable  laminated  shell 
elements  is  developed  such  that  the  high-fidelity  ANCF  tire  model  can  be  integrated  into  general 
multibody  dynamics  computer  algorithms  for  ground  mobility  simulation. 

To  this  end,  the  tire  cross  section  geometry  is  imported  from  the  tire  cut  section  and  data 
points  are  interpolated  by  a  cubic  smoothing  spline  for  extracting  the  nodal  position  and  slope 
coordinates.  As  shown  in  Fig.  4,  the  tire  cross-section  is  divided  into  the  tread,  sidewall,  and 
bead  sections.  The  number  of  layers,  cord  angles  of  layers,  material  properties  are  provided  in 
each  section  to  create  the  tire  model  data  as  shown  in  Fig.  5.  The  tread  section  consists  of  a 
carcass  ply,  two  steel  belts,  a  belt  cover,  and  tread  blocks.  The  carcass  ply  and  steel  belt  are 
modeled  as  an  orthotropic  material  with  nylon  and  steel  cords  embedded  in  rubber,  respectively. 
A  rubber  layer  is  considered  between  the  upper  and  lower  steel  belts  and  between  the  carcass  ply 
and  the  lower  steel  belt.  The  sidewall  section  is  modeled  by  two  carcass  plies  and  a  rubber  that 
lies  in  between.  The  bead  section  is  modeled  by  two  carcass  plies,  a  steel  belt,  and  a  rubber  as 
shown  in  Fig.  5.  Having  determined  the  cross-section  property,  the  three-dimensional  tire 
geometry  is  generated  by  rotating  the  tire  section  model  and  the  nodal  coordinates  of  the  ANCF 
tire  model  are  created  as  summarized  in  Fig.  4. 

The  tire  air  pressure  is  220  kPa  that  is  considered  by  the  normal  distributed  load  applied  to  the 
inner  surface  of  the  tire.  The  penalty  approach  is  used  for  modeling  the  normal  contact  force  at 
each  node  in  contact.  The  load-deflection  curve  is  important  for  characterizing  the  fundamental 
structural  properties  of  tires.  The  lateral  and  vertical  deflections  of  the  ANCF  tire  model  shown 
in  Fig.  6  are  compared  with  the  measurement  results  in  Figs.  7  and  8.  It  is  observed  from  these 
figures  that  the  local  tire  deflections  are  well  predicated  in  both  lateral  and  vertical  directions  for 
the  various  wheel  loads.  Furthermore,  the  lengths  of  the  contact  patch  in  the  longitudinal  and 
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lateral  directions  also  agree  well  with  those  of  the  measurement  results  as  shown  in  Fig.  9  and  10. 

5.  SUMMARY  AND  CONCLUSIONS 

In  this  study,  a  continuum  mechanics  based  shear  deformable  shell  element  of  the  absolute  nodal 
coordinate  formulation  (ANCF)  is  generalized  to  a  laminated  shell  element  for  application  to  the 
modeling  of  fiber-reinforced  rubber  (FRR)  structure  of  the  physics-based  ANCF  tire  model.  It  is 
shown  that  the  complex  deformation  coupling  exhibited  in  fiber-reinforced  composite  materials 
can  be  automatically  considered  in  the  shear  deformable  laminated  shell  element  using  the 
continuum  mechanics  approach.  Furthermore,  the  element  lockings  are  systematically  eliminated 
by  the  assumed  natural  strain  and  enhanced  strain  approaches,  thereby  leading  to  a  locking-free 
shear  deformation  ANCF  laminated  shell  element.  The  benchmark  problem  of  the  extension  and 
twisting  coupling  of  two-layer  laminated  plate  is  used  to  validate  the  ANCF  laminated  shell 
element.  The  numerical  results  are  in  good  agreement  with  those  predicated  by  the  analytical 
model  based  on  the  classical  lamination  theory  and  the  laminated  solid  shell  element  in  ANSYS. 
Furthermore,  using  the  ANCF  laminated  shell  element  developed,  a  physics-based  ANCF  tire 
model  is  developed  by  considering  the  detailed  tire  geometry  and  material  properties.  The  fiber- 
reinforced  rubber  is  considered  for  modeling  the  carcass  plies  and  steel  belts  using  the  multi¬ 
layered  laminated  shell  elements.  The  load-deflection  curves  predicted  by  the  physics-based 
ANCF  tire  model  are  in  good  agreement  with  the  measurement  results.  In  the  future  work,  the 
tangential  tire  force  model  based  on  the  distributed  LuGre  friction  model  is  integrated  into  the 
ANCF  tire  model  (ANCF-LuGre  tire  model)  for  use  in  the  transient  braking/traction  and 
cornering  force  evaluation  of  tires  in  the  context  of  multibody  vehicle  dynamics  simulation. 
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Figure  1.  Kinematics  of  shear  deformable  ANCF  laminated  shell  element 
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Figure  2.  Uniaxial  tensile  test  of  two-layer  laminate 
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Figure  3.  Twisting  of  two-layer  laminate  subjected  to  uniaxial  tensile  load 
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1.  Cut  section 


4.  Element  discretization 


5.  3D  tire  geometry 


6.  ANCF  tire  model  data 


Figure  4.  ANCF  tire  model  creation  procedure 
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Figure  5.  Physics-based  ANCF  tire  model  using  multi-layered  laminated  shell  element 
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Figure  6.  Deformed  shape  of  tire  cross  section 
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Figure  7.  Lateral  deflection  of  tire  for  various  wheel  loads 
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Figure  8.  Vertical  deflection  of  tire  for  various  wheel  loads 
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Figure  9.  Longitudinal  contact  patch  length  for  various  wheel  loads 


This  material  is  declared  a  work  of  the  U.S.  Government  and  is  not  subject  to 

copyright  protection  in  the  United  States.  Approved  for  public  release;  distribution  is  unlimited. 


Lateral  contact  patch  length  (mm) 


UNCLASSIFIED:  Distribution  Statement  A.  Approved  for  public  release.  #26261 


Figure  10.  Lateral  contact  patch  length  for  various  wheel  loads 


This  material  is  declared  a  work  of  the  U.S.  Government  and  is  not  subject  to 

copyright  protection  in  the  United  States.  Approved  for  public  release;  distribution  is  unlimited. 


